Critical velocity of a mobile impurity in one-dimensional quantum liquids 
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We study the notion of superfluid critical velocity in one spatial dimension. It is shown that 
for heavy impurities with mass M exceeding a critical mass Mc, the dispersion develops periodic 
metastable branches resulting in dramatic changes of dynamics in the presence of an external driving 
force. In contrast to smooth Bloch Oscillations for M < Mc, a heavy impurity climbs metastable 
branches until it reaches a branch termination point or undergoes a random tunneling event, both 
leading to an abrupt change in velocity and an energy loss. This is predicted to lead to a non-analytic 
dependence of the impurity drift velocity on small forces. 



The study of the effect of impurities constitutes a ma- 
jor theme in condensed matter physics. In the solid state 
impurities are typically regarded as immobile, acting as 
static perturbations. Prime examples include the resid- 
ual resistance of metals and the Kondo effect. In other 
(typically liquid) states mobile impurities play a vital 
role, with the influence of ■^He impurities on the trans- 
port properties of superfluid *IIe being perhaps the best 
studied example [1, 2] . 

The traditional goal has been to understand the effect 
of foreign particles or defects on an otherwise pure sys- 
tem. However, the converse problem - to determine the 
effect of the surrounding medium on the motion of impu- 
rities - is increasingly relevant to the study of ultracold 
atomic gases. There are several ways in which mobile 
impurities appear naturally in these systems: (i) a small 
fraction of atoms may be transferred to a different hy- 
perfine state, creating a dilute gas of distinguishable im- 
purities of the same mass [3-5]; (ii) mixtures of different 
atomic species {e.g., ^^Rb and ^^K in Ref. [6]) with one 
dilute component give a route to impurities of mass M 
different from the mass m of the host particles; (iii) Ions 
of Yb"^ , Ba"*" or Rb"*" have been placed in a Bose- Einstein 
condensate of neutral ^^Rb atoms [7, 8] . In each of these 
situations, impurity dynamics may be investigated with 
the help of external forces that act selectively on the im- 
purity. 

Owing to recent advances in cold atom trapping, sys- 
tems in reduced dimensionalities are now well within ex- 
perimental reach [5, 7, 8], bringing to light a wide range 
of rich and peculiar phenomena. In one spatial dimen- 
sion (ID) Refs. [9, 10] established that even in the limit of 
strong interactions, there exists a timescale beyond which 
an impurity with quadratic dispersion, E{P) = P'^ /2M* , 
may propagate in a quantum fluid. In this Letter we fo- 
cus on the full dispersion E{P) of the resulting motion, 
which possesses several remarkable features not present 
in higher dimensions. The most prominent property is its 
periodicity, which is not evident from a quadratic, small 
P expansion, and may be explained as follows. The mo- 
tion of A'^ particles on a ring with circumference L re- 
sults in the total momentum being quantized in units of 
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FIG. 1. Dispersion law for an impurity in a ID quantum 
fluid for various impurity masses, (a) M < Mc, E{P) is 
a smooth function of P and H (inset, Eq. (2)) has a single 
minimum, (b) M > Mc, there exists a metastable minimum. 
An impurity driven by a force climbs a metastable branch 
until tunneling or H loses its minimum. Each cycle releases 
energy e to the system, (c) M ^ Mc, many minima co-exist. 



2TThN/L ~ 2'Khn, where n is the one-dimensional density 
of the fluid. On the the other hand, in the thermody- 
namic limit N,L ^ oo, N/L = n, this motion costs 
no kinetic energy due to infinite total mass Nm. Thus, 
instead of paying macroscopic energy {2Trhn)^ /2M* to 
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deliver this momentum to impurity motion, the system 
channels it to the superflow of the background liquid at 
no energy cost. As a result, the dispersion of the dressed 
impurity is a periodic function E{P) = E{P + Inhn). 
Owing to the quadratic dispersion near P = 0, impurity 
motion costs less energy for a given momentum than the 
softest excitations of the liquid, i.e., phonons with lin- 
ear dispersion. This fact implies that impurity motion, 
while possessing a periodic dispersion, also defines the 
lower edge of the many-body continuum, above which 
(shaded grey. Fig. 1) there exist excitations comprised 
of the moving impurity and phonons [11-14]. Explicit 
examples demonstrating these features are furnished by 
models solvable via Bethe Ansatz [15-17]. A dramatic 
consequence of the dispersion periodicity is the possibil- 
ity of observing Bloch oscillations of a driven impurity in 
ID in the absence of a lattice [18, 19]. 

The same picture of momentum transfered to the fluid 
at no energy cost leads to the prediction of vanishing Lan- 
dau critical velocity in ID [20, 21]. It was shown that an 
infinitely massive object moving with any nonzero veloc- 
ity with respect to the background nucleates phase slips 
corresponding to one quantum 2TThn of total momentum 
transferred to the fluid. The finite rate of this process 
can then be related to the energy and momentum dis- 
sipation in the liquid, precluding superfluidity. On the 
other hand, a mobile impurity with finite mass M can 
propagate without dissipation at zero temperature with 
velocity given by slope of the dispersion V = dE/dP. 
Hence, the maximum slope of the dispersion defines a 
non-zero critical velocity for light mobile impurities. 

The main result of this work is the prediction that the 
above regimes of heavy and light impurities are separated 
by a quantum phase transition taking place at a critical 
mass Mc, given by Eq. (3). We show that above the crit- 
ical mass the many-body ground state in a total momen- 
tum sector P has cusps at P = ±hTrn, ±3/i7rn, . . ., while 
the impurity dispersion develops a swallowtail structure 
with mctastable branches shown in Fig. 1. In this case 
the true thermodynamic critical velocity Vc, given by the 
maximal slope of the many-body ground state, should be 
distinguished from the dynamic critical velocity Vc ob- 
tained from the maximal slope of the metastable branch. 
It can be shown that for M ^ the thermodynamic 
critical velocity is inversely proportional to M consistent 
with the infinite mass limit of Refs. [20, 21]. 

Metastability and two definitions of critical velocity. 
For a concrete example of how a metastable impurity 
dispersion can arise, let us take our fluid to be a weakly 
interacting Bose gas. At the mean fleld level, appropriate 
for weak interactions, the condensate wavefunction devel- 
ops a phase drop $ across the impurity located at X. In 
the limit of strong coupling between gas and impurity, 
the latter gives rise to the Josephson term 

= -nKcos* (1) 

(we have set h = 1) in the impurity energy. Here Vc 
is the dynamical critical velocity which depends on the 



details of the impurity. The phase drop inevitably creates 
a background supercurrent contributing a term n$ to the 
total momentum P. The total energy of the impurity and 
the background can then be written as 

The first term is the kinetic energy of the impurity mov- 
ing with velocity V = {P — n^)/M, where M is the total 
mass of the impurity and induced depletion cloud. As we 
shall see one can assume ~ M for a sufficiently heavy 
impurity. 

The phase drop $ represents a collective coordinate 
characterizing the state of the impurity equilibrated with 
the background [19]. Its value is determined from the re- 
quirement of the minimum of the total energy (2) leading 
to the matching P — n$ = MVc sin $ of the Josephson 
current nV^sin# to the current nV experienced by the 
impurity in its rest frame. For sufficiently large mass 
M > Mc this equation has several solutions correspond- 
ing to different arrangements of the impurity and the 
background. The total energy H{P, $) in Eq. (2) can 
therefore have minima for several values of the phase drop 
$ for a given value of the total momentum P. Plotting 
the corresponding energies results in a typical swallowtail 
structure for the dispersion, see Fig. 1. 

Following the global minimum E_{P) of the disper- 
sion one sees that it develops cusps at the crossing points 
P = ±7rn, ±37rn, . . . corresponding to a first order transi- 
tion as a function of the total momentum P. In addition 
to the true ground state E- (P) there is also a metastable 
branch E+{P) corresponding to local minimum. The 
minima corresponding to E-{P) and E+{P) are sepa- 
rated by the maximum. At the termination points Pt 
the maximum merges with the local minimum at $ = $t 
and the metastable branch ceases to exist. The critical 
mass Ale for the swallowtail catastrophe can be estimated 
from the simple model Eq. (2) as 

Mc = T^n/Vc = KmicjVc) , (3) 

where K = im/mc is the Luttinger parameter of the 
liquid depending on the speed of sound c. For a strongly 
repulsive impurity in a weakly interacting background 
Mc ^ m since K 1 and c:$> Vc, justifying a posteriori 
the assumption M. M. 

For M » Mc the termination point corresponds 
to sin"I>t = ±1 and the mass- independent velocity 
dE^/dP = ±Vc. In contrast, the thermodynamic crit- 
ical velocity, defined with the help of the stable branch 
Vc = dE-/dP, behaves as 1/M resulting in zero critical 
velocity for infinitely heavy impurity [20, 21]. 

Quantum fluctuations can smear the transition by pro- 
ducing a mechanism for the decay of the metastable 
branch [14]. The metastable branch E+{P) is well de- 
fined as long as the decay rate r(P) to the lower branch 
satisfies r(P) < E+{P) - E_{P) « 2Vc|P - 7rn|. We 
show below that r(P) vanishes as a power law near the 
crossing points, r(P) oc )P — 7rn)" for P — )• tto with the 
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exponent a given in Eq. (7). Self-consistency therefore 
requires a > 1. Hence a = 1 is the condition defining the 
transition between the two classes of dispersion curves, 
and corresponds to M = [14]. 

Before proceeding we note that the qualitative fea- 
tures discussed above in no way depend upon the specific 
model Eq. (1) and occur (for example) in the dispersion 
relation of a particle moving in a Bose gas described by 
the Gross Pitacvskii equation, valid for arbitrary cou- 
pling between the gas and impurity [22]. 

Depleton model. We now introduce a framework to 
describe the periodic impurity dispersion, the metastable 
branch, and its decay in the general case. The moving 
depleton, i.e. impurity dressed by the liquid depletion, is 
characterized by a number of particles N expelled from 
its core in addition to the phase drop The existence 
of the two collective slow coordinates is associated with 
the presence of the two conservation laws: number of 
particles and momentum [19]. The depleton's dynamics 
is then specified by the following Hamiltonian 
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{P-n^f + ljN + Hd{^,N). (4) 



Here is the ^-periodic energy function of the deple- 
tion cloud, fi is the chemical potential of the liquid in 
the absence of the impurity and A4 = M — mN. In 
the limit of strongly repulsive impurity the dynamics of 

is frozen and Hd{^, N) reduces to the simple Joseph- 
son form Eq. (1). If the mass M. is sufficiently large, the 
Hamiltonian (4) (being almost a periodic energy function 
of <&) possesses many metastable minima in addition to 
the absolute one. 

We wish to determine the tunneling rate from the 
above mentioned metastable minima to the ground state 
branch at the same momentum P. Such tunneling is ac- 
companied by a change in the phase drop A$ and num- 
ber of depleted particles AiV, despite the momentum of 
both states being equal [23]. Therefore one must know 
the Lagrangian governing the dynamics of the collective 
variables $ and N. As shown in Ref. [19], this requires 
the introduction of, and coupling to, the phonon subsys- 
tem. The latter may be described by small deviations of 
the density, p{x,t), and velocity, u{x^t), fields from their 
unperturbed values. It is convenient to introduce the su- 
perfluid phase (p{x,t) and the displacement field ^}{x,t) 
defined by zt = dx^p/m and p = dx^/n. 

The effective action governing these harmonic degrees 
of freedom is that of the Luttinger liquid [25, 26] 
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while their coupling to the depleton degrees of freedom 
takes a universal form [19, 24], 
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d{X,t)/TT - N^p{X, t) 



(6) 



It is the phonon subsystem which must supply the macro- 
scopic momentum nAcl? and take away AA'^ particles from 



the depleton. Near the point P = irn, the phonons have 
only a small amount of energy e with which to do this. 
This low-probability event, or 'under-the-barrier' evolu- 
tion of the phononic system, can be described in imagi- 
nary time by performing the Wick rotation t — >■ — ir in 
Eqs. (5), (6). We have delegated details of the tunnel- 
ing calculation to a short section in the Supplemental 
Material [22]. We find for the tunneling rate 
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Writing e = 2Vc \P — TTn\ gives the momentum depen- 
dent tunneling rate advertised earlier. 

The logarithmic behavior of the tunneling action [22] 
is valid only on the long time-scale t > //c, implying the 
validity of Eq. (7) is restricted to the regime of small 
£ close to P = nn. In such a case the velocities of 
the upper and lower branches of the dispersion are es- 
sentially opposite and close to ±Vc, while the change in 
number of particles AN, being an even function of veloc- 
ity, is negligible. As a result one arrives at the exponent 
a = 2if (A$/27r)^ — 1. This single collective coordinate 
limit, i.e. AN = 0, agrees with the previous result of 
Ref. [28] (cf. Eq. (7.20) therein) based on a phenomeno- 
logical model of tunneling between two minima in the 
presence of coupling to an environment. Their results 
also show that a finite temperature T acts as a low en- 
ergy cut-off, namely F ~ T", e <^ T (see Eqs. (5.30), 
(7.18) of Ref. [28]). 

Higher order processes involve tunneling back and 
forth between the two minima that become degenerate 
for e = 0, leading to an effective spin-boson or anisotropic 
Kondo description [14]. The inclusion of such processes 
is only necessary if a < 1, meaning that tunneling is a 
relevant perturbation. This leads to the rounding out of 
the cusp in the dispersion (P) at P = nn, and the dis- 
appearance of the metastable branch. Thus a = 1 is the 
exact condition determining the critical mass M = Mc. 

In the infinite mass limit with fixcid velocity V = 
P/M Eq. (4) becomes essentially a washboard potential. 
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adjacent minima goes to zero in this limit). Thus the 
change A$ between the adjacent minima is exactly 27r, 
while the energy difference is e = 2'imV. As a result, 
one finds a = 2K — 1 and the power law dependence of 
the tunneling rate F ~ I/^at-i |'29] on the velocity in full 
agreement with the result of Ref. [20, 21]. 

In the limit of large energy detuning, i.e., away from 
P = nn, the tunneling rate does not take the form (7). In 
this region the barrier is small and the impurity emerges 
from under it with only slightly less energy, despite the 
large detuning between the branches (see Fig. 1 insets). 
Therefore most of the descent from the upper branch 
corresponds to viscous relaxation in real time and not an 
under-the-barrier tunneling trajectory. Thus, near the 
termination momentum Pt ~ MV^, the potential may 
be approximated by a cubic polynomial. The dissipa- 
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FIG. 2. Impurity velocity (solid) in units of the dynamic 
critical velocity K as a function of time in units of the period 

T = 2Tin/F. The drift velocity is denoted by Vb (dashed) and 
is displayed in the inset as a function of the applied force. Due 
to quantum tunneling near P — nn, the drift velocity goes to 
zero as a power law with exponent l/(a + 1) (solid, inset) 
given by Eq. (9) for F < Fc and approaches the semi-classical 
result for F > Fc (dashed, inset) [22]. 



tive tunneling through such a potential was described in 
Ref. [30], leading to the escape rate F cx , where 



S = 2K 
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(8) 



We see that for both large and small energy detuning, 
Eqs. (7), (8) the tunneling action scales with the Lut- 
tinger parameter K. 

One way to explore the metastable branches of the dis- 
persion is by applying an external force F to the impu- 
rity and studying the ensuing dynamics. For sufficiently 
strong force the driven impurity overshoots the ground 
state branch and follows the metastable branch until it 
either tunnels or reaches the termination point Pt (see 
Fig. 1). The energy dissipated per cycle, e, must be sup- 
plied by external force: FVd = s/T resulting in the drift 
velocity Vb = e/2-Kn. Here T = 2'Kn/F is the period of 
the motion. 

Tunneling is negligible when FT <C 1, which occurs 
if the force is sufficiently large. In such a case the ex- 
ternal force drives the impurity all the way up to the 
termination point Pf It then loses the energy minimum 
and relaxes to the lower branch, emitting phonons. As 
a result, the energy dissipated per cycle and thus drift 
velocity are essentially independent of the applied force. 
This is in stark contrast to the case of M < Mc, where 
the phonon emission due to dipole radiation gives rise to 
a linear in F drift velocity and hence a finite impurity 
mobility ct = Vd/P [19]. 

The impurity velocity takes the form of a saw-tooth 
trajectory, experiencing a slow build-up and a sudden 
drop when the velocity reaches 14, see Fig. 2. As a 
function of the mass, Vt> saturates to the critical ve- 



locity in the limit M ^ co. This is most easily seen 
by noting that in such a limit the energy drop occurs 
between two adjacent and essentially parallel parabolas: 
e = 2TmVT, = PV2M - (P - 2-Knf I2M 2-KnV^. On 
the same grounds, the amplitude of the saw-tooth oscil- 
lations is a decreasing function of M going to zero in the 
infinite mass limit. 

Quantum effects become dominant only if the time 
spent on a metastable branch is comparable with the 
inverse tunneling rate: FT > 1. Thus, the effect of quan- 
tum tunneling plays a role only for very small forces. In 
this limit, the impurity is likely to tunnel at a momen- 
tum immediately above P = 7rn, right when it enters the 
metastable branch. Using the rate Eq. (7) to calculate 
the tunneling- averaged energy drop (e), we find for the 
drift velocity [22] 
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The non-analytic dependence of the drift on small forces 
is the result of the cusp of the ground state energy at 
P = Tin and implies the divergence of the mobility cr = 
lii/F as P — across the M = transition (see inset 
of Fig. 2). At finite temperature T the mobility is not 
given by (7 ^ -y--a g^g niight be expected due to the finite 
tunneling rate F ~ T" near the cusp. Strictly speaking, 
it is given by a ~ due to two-phonon thermal Raman 
scattering [18, 19, 31] which results in velocity saturation 
below Vc as P -T- 0. In practice, therefore, it is necessary 
to work at sufficiently low temperatures and finite drives 
to overcome such saturation and allow the impurity to 
explore the full range of momentum where metastable 
branches become important. 

In conclusion, we have shown that the dispersion re- 
lation of a mobile impurity undergoes a sudden change 
when the mass of the impurity exceeds a certain critical 
vahie. The latter depends on the interaction parameter 
within the host liquid and thus may be varied experimen- 
tally. The predicted transition can be probed by applying 
an external force and monitoring either impurity velocity 
or excitations of the host liquid. The change of disper- 
sion also affects the oscillation frequency and dissipation 
of the impurity in the trapping potential. With investi- 
gations of impurity motion in ID atomic gases now un- 
derway [5, 7, 8], we hope that these predictions will guide 
future experiments. 
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I. SUPPLEMENTAL MATERIAL 

Impurity coupled to a weakly interacting liquid - To 
model the heavy impurity in the presence of a weakly 
interacting superfluid at T = 0, we employ the Gross- 
Pitaevskii equation (GPE) using a point description of 
the impurity wavefunction with coordinate X, 
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+ G6{x- X)]"^. (10) 



We look for a solution of the traveling wave form 
'b{x,t) = - Vt) and X{t) = Vt. The usual way of 
dealing with a delta function perturbation is to attempt a 
solution constructed from two homogeneous (i.e., G = 0) 
solutions that must satisfy the proper boundary condi- 
tions ai X — X. This strategy is facilitated by the fact 
that the bare GPE (G = 0) admits a one-parameter fam- 
ily of solutions 
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known as grey solitons [32, 33]. They can be visualized 
as a density dip moving with velocity V , having a core 
size I = ^(1 — F^/c-^)^^/^, where c = ^Jng/m is the speed 
of sound, n = 11/ g is the ID density and ^ = 1/mc is the 
healing length (the condition of weak inter-particle cou- 
pling implies ^ 1). Thus, by appropriately matching 
two solitonic solutions at x = Vt, one formally solves 
Eq. (10) with 'I'(a; - Vt) given by 

where y = x — Vt. Here we introduced two "matching pa- 
rameters" Xq and $0 which must be chosen to satisfy the 
two boundary conditions: ^-(0+) = *(0"), ^'I^Zq- = 
2mG'^{0). This results in the following two equations 
for $0 and z — tanh(a;o/0 
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where we have parameterized the velocity V by the soli- 
tonic phase $s via V/c = cos^. 

Equations (13), (14) permit a solution only for V < Vc 
where Vc = Vc{G) < c is the dynamic critical velocity 
that depends only upon the parameter G [34, 35]. This 
may be seen by plotting the right and left hand sides of 
Eq. (14). While the left hand side is bounded by a max- 
imum, the right hand side grows quadratically with z 
and therefore the solution exists only for a limited range 
of $s- For this reason we choose to parameterize the 
solution, Eq.(12), by the total phase drop across the im- 
purity, $ = $s — which permits a solution for any $. 
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FIG. 3. Superfluid density (top panel) and phase (bottom 
panel) profiles in the presence of an impurity propagating 
with velocity V/c = 0.66 and coupling G/c — 0.25 (solid 
lines). The repulsive impurity tends to deplete the host liquid, 
pushing away N particles and inducing a phase drop $ across 
it. The dashed lines correspond to the double soliton solution 
Eq.(ll) used to construct *(j; — Vt). 



Upon solving Eqs. (13), (14) one thus finds the periodic 
relations z — z{^,G/c) and $s = 4>s($,G/c) (i.e., they 
are invariant with respect to $ — >■ 4> -I- 27rj for integer j). 

We may use these relations to compute the dressed 
impurity energy E — ^AIV^ + and momentum P = 
MV + P(i, where E^ and are the corresponding deple- 
tion cloud contributions given by 



Ed ~ I dx 
4 
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(15) 



(16) 



Here we introduced the number of expelled particles 
~ mc ^(1 ~ •^)- Written in this way, the equilib- 
rium expressions E{^ + 27rj) = E{^) and F($ -I- 27rj) — 
2Trnj+P{^) are functions of the total phase drop $ across 
the impurity, whose functional forms depend only on the 
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FIG. 4. Top Panel; Variation of the impurity dispersion 
across the M = Mc transition, exhibiting the swallowtail 
catastrophe. Here Vc/c = 0.1, K = lOvr and Mc = 12n/c 
(see Eq. (18)). Bottom panel: Momentum P(4') for different 
masses, see Eq. (17). For M > Mc there persists two extrema 
which merge and disappear at M = Mc- 
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FIG. 5. Top panel: Impurity dispersion at weak coupling with 
parameters Mjra = 60, K = IOtt, G/c = 0.5 (green), G/c = 
0.2 (blue) and G/c = 0.05 (red). Bottom panel: Momen- 
tum P as a function of $ having the same parameters listed 
above. The G — >■ limit is well approximated by the gray 
soliton configuration (dashed): P — n($ — sin<I>) + Afccos*, 
except near $ = lixj for integer j, where the momentum is 
necessarily given by P = 2-Knj. 



parameter G/c. Alternatively, one may solve for $ = 
<i>(P) with the help of Eq. (16). Substituting the resulting 
relation into the energy function, i?(<I>(P)), one obtains 
the periodic dispersion relation E[P~) — E{P + 27rn). 

The above point description of the impurity is mean- 
ingful only if its effective mass is sufficiently large. This 
can be achieved either by a large bare impurity mass or 
strong coupling to the liquid. The latter tends to localize 
the impurity near the minimum of the self-induced deple- 
tion cloud density, see Fig. 3, while the former leads to 
a small de Broglie wavelength. To allow ourselves the 
generality of discussing a wide range of impurity masses, 
and for simplicity of illustration, we momentarily focus 
on the regime of strong coupling G ^ c. 

Strong impurity coupling - We determine the depen- 
dence of z and $s on $ to leading order in c/G. From 
Eqs. (13), (14) we find z{<P) ~ ^cos^f and $«($) ~ 



TT — ^ sin $ and hence, 

P = n<i>+ (^M - Vcsm<t>. (17) 

Here the critical velocity is given by Vc = ^ ^ c [34, 35] , 
and the $ dependence of the velocity acquires the Joseph- 
son form: V{^) = VcSin$. The appearance of the critical 
mass occurs when P = P{^) is satisfied by multiple val- 
ues of <i>(P). The functional form of P($) for M > Mc 
takes an "J\f" shape, implying the existence of two ex- 
trema in the range $ S (0, 2tt), see Fig. 4. As M crosses 
Mc from above, these extrema merge and disappear at 
P = n$ = nn. This observation allows to derive the 
critical mass within the above framework, 

dP 2n / G\ mK c , , 
— = =^ Afc = — 1 + - « , (18) 

a$ <I.=^ C \ C J IT Vc ^ ' 

in agreement with the parametric dependence Eq. (3) of 
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the main text. 

Weakly coupled, heavy impurity - The case of weak im- 
purity coupUng does not lend itself to convenient analyti- 
cal analysis. For given G/c < 1, we thus solve the bound- 
ary Eqs. (13), (14) numerically to produce the disper- 
sion law E{P), see Fig. 5. For comparison to the strong 
coupling case, we also plot P(i>) for various couplings 
G/c < 1. In the limit G ^ 0, the impurity is dressed 
by an essentially unperturbed gray soliton {i.e., the two 
solitons of Fig. 3 are separated by a very small distance) . 
This approximation breaks down near $ = 27rj for in- 
teger j since we must have P = 2TTnj. Consequently, 
there is a sharp deviation from the soliton configuration 
(dashed line in Fig. 5) near these points, making the func- 
tional form of P(^>) difficult to characterize analytically. 
As a result of the numerical analysis, we see that even in 
the regime of weak coupling the appearance of multiple 
momenta extrema give rise to the observed swallowtail 
catastrophe of the energy dispersion and is thus not re- 
stricted to the Josephson limit. 

Quantum tunneling and non-analytic drift - We wish 
to calculate the tunneling-averaged energy drop (e) as 
a function of the (small) force F = P/t using the mo- 
mentum dependent tunneling rate r(P). To this end we 
establish the validity of the result (7) in the main text. 

The tunneling rate r(f) is given by the exponentiated 
action <S'(f ) evaluated on the trajectory where the system 
tunnels for an imaginary time f from one minimum to the 
other. For a > 1 where tunneling is an irrelevant per- 
turbation, we may safely neglect higher order processes 
involving tunneling back and forth between the minima 
and focus on the single tunneling event. Starting from 
Eqs. (5), (6) of the main text (after the Wick rotation 
t — > —ir) one integrates the imaginary time Luttinger 
liquid action with linear coupling to arrive at the tunnel- 
ing action 

S{f) = -^ j dxdx'dTdT'jl^^{f)D{x-x',T-T')J^,^r'{f) 

(19) 

where 4,^(f) = 5{x-X{t)){^,N) = 5{x- X{T))[d{T) - 
(5(t — f)](A$, AA^) describes the tunneling trajectory be- 
tween minima having configurations which differ in phase 
by A$ and number of expelled particles AA'' and D{x, r) 
is the matrix imaginary time phononic propagator. The 
latter is determined by inverting the matrix operator 
of the quadratic Luttinger liquid action, Eq. (5) of the 
main text, in the form Sll = — ^ J dxdrx^ D~^Xj where 



Substituting the above tunneling trajectories into 
Eq. (19) with the imaginary time propagator 



D{x,t) 



1 

47r 




In I 



—In I 



cT-\-ix 



(20) 

we find S{t) = (1 -I- a)h\{cT/l) where / > [mc) ^ 
is a short distance cut-off beyond which the point- 
particle description of the depleton breaks down and 



a = 2K 



( ^\ ( 

\ 2tt ) ^ \ 2K I 



1. The X dependence of the 

action is neglected because it enters as X{t) — X{t') ~ 
V{t -t') <^c{t -t'). 

The energy-dependent tunneling rate is found with 
the help of the Fourier transformation, r(£) ~ 
/ die~'^(**)e"^*, where the analytical continuation to real 
energy ie — >■ e is accompanied by the Wick rotation r 
it. Substituting S{t) given above and using the identity 

/ die*^*(it)-(i+") = ^^^£", we arrive at Eq. (7) of the 
main text. 

To compute (e) one must know the probability to tun- 



nel per unit time: Pit) = r(t)Exp - dt'V{t') , which 

follows from the assumption that the process is Marko- 
vian. By parameterizing the average over time by an 
average over momenta one finds 



(£)= j ' dPe{P)^-^Exp - I dP' 

= / dxe{P{x))e-='. 
Jo 



MP') 



(21) 
(22) 



Here we introduced the dimensionless variable x{P) = 
J^^dP'T{P')/F along with = dP'T{P'). Due to 
the exponential decay of the kernel in Eq. (21), the in- 
tegral is cut at Xc = min{l, Ft^/F}. As a result, one 
obtains the approximation (e) ~ XcS{Pc), where Pc sat- 
isfies Xc = dPr{P)/F. For F < we have Xc = 1 
and hence Pc may be solved by writing Pc = mr + SP for 
small, yet undetermined, SP. Expanding F near mr as 
r = r^[(P - n7r)Vc/mc2]«, we find 5P = ^{F/Fc)^ 

and Fc = Jt.t!. w for large M. Thus, 



Vc(l-Ka) 



{s) = mc^ 



establishing Eq. (9) of the main text. 
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